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Abstract 

We study the effect of the cosmological constant on the statis- 
tical properties of gravitational lenses in flat cosmologies (Qq + 
Ao = 1). It is shown that some of the lens observables are 
strongly affected by the cosmological constant, especially in a 
low-density universe, and its existence might be inferred by a 
statistical study of the lenses. In particular, the optical depth of 
the lens distribution may be used best for this purpose without 
depending much on the lens model. We calculate the optical 
depth (probabilty of a beam encountering with a lens event) for 
a source in a new picture of galaxy evolution based on number 
evolution in addition to pure luminosity evolution. It seems that 
present day galaxies result from the merging of a large number of 
building blocks. We have tried to put limit on the cosmological 
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constant in this new picture of galaxy evolution. This evolu- 
tionary model of galaxies permits larger value of cosmological 
constant. 
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1 Introduction 

In the early history of modern cosmology the cosmological con- 
stant A was invoked twice. First by Einstein to obtain static 
models of the universe. Next by Bondi, Gold and Hoyle to re- 
solve an age crisis and to construct the universe that satisfied 
the " Perfect Cosmological Principle", i.e., the hypothesis that 
the universe appears the same at all times and places. In both 
instances the motivating crisis passed and the cosmological con- 
stant was put aside. However the cosmolgical constant remains 
a focal point of cosmology and of particle theory. The former 
because today an understanding of a wide range of observations 
seem to call for a cosmological constant. The latter because 
in the context of quantum-field theory a cosmological constant 
corresponds to the energy density associated with the vacuum 
and no known principle demands that it vanish. 

Various aspects of recent observations suggest reconsidera- 
tion of a nonvanishing cosmological constant (Krauss &; Turner 
1995). These includes the age of universe once again, the forma- 
tion of large-scale structure (galaxies, clusters of galaxies, super- 
clusters) and the matter content of the universe as constrained 
by dynamical estimates, Big Bang Nucleosynthesis and X-Ray 
observations of clusters of galaxies. 

One of the possibilities of detecting A is the gravitational 
lensing technique. To use the gravitational lensing as a tool for 
the determination of cosmological parameters either by a de- 
tailed study of specific lens systems or through statistical anal- 
ysis of a sample of lenses has been frequently discussed (Refsdal 
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1964; Press & Gunn 1973). It has been pointed out that the 
expected frequency of multiple imaging lensing events is quite 
sensitive to cosmological constant (Fukugita, Futamse & Kasai 
1990; Turner 1990). To put a limit on A, we need to first cal- 
culate the expected number of multiple image gravitational lens 
systems (produced by the known galaxy population) to be ex- 
pected in a particular quasar sample with a known distribution 
of redshifts. This then has to be compared with the observed 
frequency of lens systems found. 

We have used a unifying model of galaxy formation which can 
answer two questions raised by the recent data on high - redshift 
galaxies. These concern the ages of, and star formation history 
in the distant radiogalaxies, and the nature of the large number 
of field galaxies revealed by faint galaxy counts. This new model 
is based on the strong number evolution in addition to pure 
luminosity evolution of the galaxies (Volmerange & Guiderdoni 
1990). 

In Section 2, we write down the New Luminosity Function 
(NLF) which has strong number evolution in addition to pure 
luminosity evolution. In Section 3, we present a new calculation 
of the total multiple image lensing cross-sections for a galaxy 
in the Singular Isothermal Sphere (SIS) approximation (Turner, 
Ostriker & Gott 1984) (hereafter TOG) based on the new galaxy 
luminosity function, velocity-luminosity correlation and veloc- 
ity dispersion. In Section 4 we write down the basic equations 
for the statistical properties of lenses for each of the galaxy mod- 
els. The quasar luminosity function for the BSP sample (Boyle, 
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Shanks & Peterson 1988) is described in Section 5. Section 6 
contains the number of lensed quasars in the BSP sample for 
the Press-Schecter Luminosity Function (PSLF) and the New 
Luminosity Function (NLF). In Section 7 we discuss the results. 

2 New Luminosity Function 

In 1990, Volmerange and Guiderdoni, proposed a unifying model 
to explain faint galaxy counts as well as observational properties 
of distant radiogalaxies. This new model of galaxy evolution 
is based on number evolution in addition to pure luminosity 
evolution. Present day galaxies result from the merging of a 
large number of building blocks and the comoving number of 
these building blocks evolves as (1 + z) 1 - 5 . It is argued that the 
present luminosity function is the well known Press - Schecter 
Luminosity Function (PSLF) 

$(L, z = 0) = f(L/L*) a exp(-L/L*)dL/L, 

with being the characterstic luminosity at the knee and (jf 
a characterstic density. These values are fixed in order to fit 
the current luminosities and densities of galaxies. Then at high 
z, the comoving number density follows the New Luminosity 
Function (NLF) 

$(L, z)dL = (1 + z) 2,1 <S>(L(l + z)\ 0)dL 

It is seen that the value 77 = 1.5 gives a fair fit to the data on 
high redshift galaxies. The functional form has the following 
properties: 
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(i) Self-similarity as suggested by the classical Press- Schecter 
(1974) prescription. 

(ii) conservation of the total comoving mass density. 

(iii) evolution of the comoving number density as (1 + z) v and 
of the knee of the function as L*(z) = L*(0)(1 + z)~ v 

3 Lensing Cross-Section For SIS Galaxies 

The probability of a beam encountering a lensing object is gov- 
erned by the parameter F (TOG 1984) which represents the 
effectiveness of distributed cosmic matter in producing double 
images. The parameter F is given by 

where no is the present comoving number density of lensing 
galaxies and v is the line of sight velocity dispersion of mat- 
ter in the galaxy. The Hubble constant Hq is normalized by Hq 
= 100 h km s _1 Mpc _1 , where 0.4 < h < 1.0. In this section we 
evaluate F from the statistics of local galaxies, assuming that 
their properties persist out to distant galaxies. The relationship 
between v and the luminosity is known empirically, i.e., L oc v 4 
for elliptical galaxies and L oc v 2 ' 6 for spiral galaxies in the B 
band. Using the luminosity function given by Volmerange and 
Guiderdoni (1990) 
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${L,z)dL = (l + z) 2r V* 



T ia 
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n L (0) = / <S>(L,z)dL 
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Using eq.(l), eq.(2), eq.(3), eq.(4) we get 
16tt 3 / 4 



F = 



which may be written as 



viT (a + i + l) (1 + zf 1 '^ 



F = F*{l + zf {1 -^ 



where 



F* = 



16tt 3 



favjT (a + - + 1 
V 7 



(2) 
(3) 

(4) 



(5) 



(6) 



(7) 



here v* is the velocity dispersion corresponding to the character- 
istic luminosity L* and 7 = 4 for elliptical galaxies and 7 = 2.6 
for spiral galaxies. According to Fukugita and Turner (1991), 



(j), = (1.56 ±0.4) x lO- 2 h 3 Mpc- 3 
a = -1.1 ± 0.1 



(8) 
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The values of the parameter F* adopted by them are summa- 
rized in Table 1. 

It has been argued that the existence of a \f3/2 correction 
factor for the dark matter velocities disperion used by Fukugita 
and Turner is incorrect (Fukugita & Turner 1991). Better dy- 
namical models (Kochanek 1993; Kochanek 1994; Brimer & 
Sanders 1993) show that the assumption of ^3/2 is incorrect, 
and is supported neither by galactic dynamics nor by the ob- 
served separations of gravitational lensing. The best fit esti- 
mate is v* = 220 ±20 km/sec for E + SO galaxies, and it is well 
constrained because the average image separation is a strong 
function of velocity dispersion (Kochanek 1996). The new val- 
ues of F* obtained by using the parameters given by Kochanek 
(Kochanek 1996) are summarized in Table 2. 

0, = (1.40 ± 0.17)/i 3 x 10- 2 Mp C - 3 



a = -1.00 ±0.15 

We note that the total value of F* determined here (F* = 
0.028) is four times smaller then the value used by TOG. Naka- 
mura and Suto (1996) have given the new values of (7,^*) = 
(3.3, 175 kms^ 1 ) for E and SO galaxies and (7, v*) = (2.9, 126 kms^ 1 ) 
for S galaxies. The morphological composition is 4>*e ± 4>*so = 
0.440* and (J)*s = 0.560* where 0* = 0.26h?Mpc~ 3 and a = 
— 1.09. The values of F* wth Nakamura and Suto parameters 
are summarized in Table 3. 
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4 Basic Equations For Statistical Lensing 

To discuss the statistical properties of gravitational lenses we 
assume that the universe is well approximated by the Friedmann 
- Lemaitre - Robertson - Walker (FLRW) geometry on large 
scale. We write the basic equations for the statistical lensing for 
two models of the lensing object; i.e. , (1) Point Masses, which 
are appropriate models for stellar mini - lensing or concentrated 
sources such as black holes, and (2) Singular Isothermal Sphere 
(SIS), which would model the matter distribution of an isolated 
galaxy. We use the following notation for distances: 

D L = d(0, z L ), D LS = d{z L , z 8 ), D s = d{0, z s ), 

where d(zi, z<i) is the angular diameter distance between the red- 
shift z\ and Z2 , and the arguments zl and zs are the redshifts of 
the lens and source respectively. The formulation and notation 
of TOG is basically followed in this paper. For Q a + A Q = 1, 



where A D = A/3H% and A = 87rGp vac and R Q is the present scale 
factor. 



This model can be considered to be a good approximation 
for many celestial bodies like "Jupiter", stars, black holes and 
even galaxies, when the light rays from background sources pass 





4.1. Point Masses 
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outside the deflectors. We first define the length a cr which char- 
acterizes the critical radial size of the lens such that 



2 AGMDolDls , im 

a cr = 2 n , ( 10 ) 

c z Dos 



where M is mass of the lensing object (TOG). Then the cross - 
section o for strong lensing events as defined by TOG is given 
by 

o = 7ra 2 cr , (11) 

The differential probabilty dr of a beam encountering a lens in 
traversing the path of dzL is given by 

cdt 

dr = n L (z)a—dz L , (12) 

^-f^MwliSh- (13) 

where ni(z) = nj,(0)(l + zl) 3 and n^(0) is the comoving num- 
ber density and Ql(0) = 87rGMnL(0)/^H 2 is the lens density 
parameter which is the ratio of the local lens density to the crit- 
ical density. The quantity cdt/dzi is calculated in the FLRW 
geometry to be 



cdt R 1 

dzL ~ (1 + z L ) y/n o (l + z L f + (1 - Q - A )(l + z L f ~o 

(14) 

i? is the Hubble distance (R a = c/H Q ) , £ stands for the look 
back time and Q D is the total mass density of the universe. By 
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integrating the differential probabilty along the line of sight to 
the source , we obtain the total probabilty 



T(Z ' ] = i" tl^ (15) 



4. 2. Singular Isothermal Spheres 

The singular isothermal sphere (SIS) provides us with a rea- 
sonable approximation to account for the lensing properties of a 
real galaxy. The lens model is characterized by the one dimen- 
sional velosity dispersion v. The deflection angle for all impact 
parameters is given by a = 47ru 2 /c 2 . The lens produces two im- 
ages if the angular position of the source is less than the critical 
angle /3 cr , analogous to the critical radius in the previous sub- 
section, which is the deflection of a beam passing at any radius 
through a SIS: 

(3 cr = aD LS /D OS i (16) 

Then the critical impact parameter is defined by a cr = Dol(3, 
and the cross- section is given by 



cr 



— ^<) 4 (^) 2 ' <"> 

The differential probability of a lensing event is given by using 
eqs.(2), (5), (12), (17), 

cdt 

dr = n L {0){l + z L f(j—dz L 
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dz L 



dr = 




(18) 



It turns out that for 7 = 4 the differential probabilty is the 
same for both PSLF and NLF. The total probabilty obtained 
by integrating the differential probability along the the line of 
sight to the source as in eq.(15), is plotted in Fig.l. with both 
PSLF and NLF for Q = 0.1 and A D = 0.9 . 

5 The Quasar Luminosity Function 

Boyle, Shanks & Peterson (1988) (BSP) have published a sample 
of 420 faint ultraviolet - excess quasars that extends from z = 0.3 
to z = 2.2 and is complete for magnitude brighter than 20.9 mag. 
Using this sample, they proposed a quasar luminosity function 
of the form 



where Mb(z), the magnitude at which there is a turnover in the 
power law slope, varies as a function of quasar redshift : 




-1 



lQ0A[M B -M B (z)}(a+l) _|_ lQ0A[M B -M B (z)}{(3+l) 
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M B (z) = Mg - 2.5 k L Log(l + z) 



(20) 



Here M| and ki, as well as a and (3 from eq.(19) are constant 
given by the best fit of the parameters to the data and is the 
normalizing factor. We have used model B from BSP, a pure 
power law evolution that fits the data very well. For this model, 
the values of the constants are M* B = —22.42, ki = 3.15, a = 
-3.79, p = -1.44 . 

6 The Number Of Lensed Quasars 

Besides depending upon the number density of galaxies and their 
properties, the number of quasars gravitationally lensed also de- 
pends upon the number of unlensed quasars. This is because of 
the amplification bias i.e., the brightening of lensed quasars, the 
number of quasars a few magnitudes fainter than the detection 
limit is important in predicting the number of lensed quasars 
in a flux limited sample. Fukugita & Turner (1991) calculated 
the number of lensed quasars expected for A Q cosmology taking 
into the account the amplification bias. The number-magnitude 
counts of the lensed quasars are given by the relation 



where P(A)dA is the probabilty that lensing will cause the mag- 
nitude of the quasar (images) to decrease by A; i,e., the images 
are brighter by a factor A; A = 2.5 Log A . 



roo 

N LQ (m) = rj o N Q (m + A)P(A)dA 



(21) 
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roo rz2 av fM2 

N LQ = / Q P{A)dA j zi dzr(z) — j Mi <$>{M B ,z)dM B , (22) 

where P(A) dA = 7.37 lCT - 8A dA and dV/dz is the comovmg 
volume which is given by 

dV __ 47Td 2 L c 

~dz~ ~ H (l + z) 2 y/n o (l + zf + A - Wo + Ao - 1)(1 +W 

(23) 

where di, is the luminosity distance. The A^q calculated in 
BSP sample by using PSLF and NLF for galaxy distribution 
are summarized in Table 4. The result is plotted in Fig. 2 

7 Conclusion 

An upper bound to A, results from the comparison of the statis- 
tics of quasars observed to be gravitationally lensed by inter- 
vening galaxies, with the predictions of flat cosmological mod- 
els with a nonzero cosmological constant. A flat cosmology 
with cosmological constant tends to produce more gravitation- 
ally lensed quasars then does such a cosmology with zero A be- 
cause a large A increases the volume per unit redshift of the 
universe at high redshift. This means that the relative number 
of lensed quasars for a fixed comoving number density of galax- 
ies increases rapidly with increasing A. Turning this around it 
is possible to use the observed probabilty of lensing to constrain 
A. This method has been used by various authors (Turner 1990, 
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Fukugita & Turner 1991, Fukugita et al. 1992) and recently 
Kochanek has put limit on A i.e., A < 0.65 (Kochanek 1996). 

However, limited observational data and the possibility that 
evolution effects on the population of lensing galaxies have not 
been properly taken into account suggest that this limit is still 
uncertain. We have tried to put limit on A by taking evolution- 
ary effect of galaxies. The important result can be described as 
follows: The probabilty of a beam encountering with lens event 
has been reduced in an evolutionary picture (NLF) and as a 
consequence of it, the evolutionary model permits higher value 
of A. This can be explained as follows, since there is no observed 
lensed quasars in BSP sample. Taking a typical uncertainity of 
« ±1 in the number of lensed quasars, PSLF put limit on A 
i.e., A < 0.6 with F* = 0.021 whereas NLF gives A < 0.7 with 
F* = 0.021. And NLF also gives A < 0.47 with F* = 0.028. 
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Figure Captions 

Fig. 1 z s Vs t/F* 
Fig. 2 A Vs N LQ 



Table 1. SIS velocity dispersion and F* val- 
ues (using parameters of Fukugita and Turner 
(1991)). 



Type 


Composition 


v*{km/ sec) 


p* 


E 


12 ± 2% 


(225)i| x ^3/2 


0.019 ±0.008 


SO 


19 ± 4% 


(206)1^ x ^3/2 


0.021 ±0.009 


S 


69 ± 4% 


144l? 3 ± 10 


0.007 ±0.003 


ALL 


100% 




0.047 ±0.019 
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Table 2. SIS velocity dispersion and F* values 
(using parameters of Kochanek (1996)). 



Type 


Composition 


v*{km/ sec) 


p* 


E+ SO 


43% 


220 ± 20 


0.023 


S 


57% 


144 ± 10 


0.005 


ALL 


100% 




0.028 
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Table 3. SIS velocity dispersion and F* val- 
ues (using parameters of Nakamura and Suto 
(1996)). 



Type 


Composition 


v*{km/ sec) 


p* 


E+ SO 


44% 


175 


0.016 


S 


56% 


126 


0.005 


ALL 


100% 




0.021 
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Table 4. Number of lensed quasars Nlq in BSP 
sample. 



+ A 


Nlq 


with NLF 


Nlq with NLF 


Nlq with PSLF 




p* 


= 0.028 


F* = 0.021 


F* = 0.021 


1.0 + 0.0 




0.58 


0.38 


0.43 


0.9 + 0.1 




0.65 


0.43 


0.49 


0.8 + 0.2 




0.72 


0.47 


0.54 


0.7 + 0.3 




0.81 


0.53 


0.61 


0.6 + 0.4 




0.92 


0.60 


0.69 


0.5 + 0.5 




1.06 


0.69 


0.79 


0.4 + 0.6 




1.29 


0.83 


0.97 


0.3 + 0.7 




1.59 


1.02 


1.19 


0.2 + 0.8 




2.11 


1.35 


1.58 


0.1 + 0.9 




3.13 


1.98 


2.35 


0.05 + 0.95 




4.13 


2.59 


3.09 


0.01 + 0.99 




5.83 


3.61 


4.37 


0.0 + 1.0 




6.51 


4.02 


4.88 
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Fig. 1 
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